Theoretical Model for Kramers-Moyal's Description of Turbulence Cascade 
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We derive the Kramers- Moyal equation for the conditional probability density of velocity increments 
from the theoretical model recently proposed by V.Yakhot jjj in the limit of high Reynolds number 
limit. We show that the higher order (n > 3 ) Kramers - Moyal coefficients tends to zero and 
the velocity increments are evolved by the Fokker- Planck operator. Our result is compatible with 
the phenomenological descriptions by R.Friedrich and J.Peinke M, developed for explaining the 
experiments recently done by J. Pienke. et al. ||. 
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The problem of scaling behavior of longitudinal veloc- 
ity difference U — \u(x\) — u(x2)\ in turbulence and the 
probability density function of U i.e P(U), attracts a 
great deal of attention [4-10]. Statistical theory of tur- 
bulence has been put forward by Kolmogorov [11], and 
further developed by others [12-15]. The approach is 
to model turbulence using stochastic partial differential 
equations. Kolmogorov conjectured that the scaling ex- 
ponents are universal, independent of the statistics of 
large-scale fluctuations and the mechanism of the viscous 
damping, when the Reynolds number is sufficiently large. 
However, recently it has been found that there is a rela- 
tion between the probability distribution function (PDF) 
of velocity and those of the external force (see [1] for more 
detail). In this direction, Polyakov [5] has recently offered 
a field theoretic method to derive the probability distri- 
bution or density of states in (l+l)-dimcnsions in the 
problem of randomly driven Burgers equation [16-17]. In 
one dimension, turbulence without pressure is described 
by Burgers equation (see also [18] concerning the rela- 
tion between Burgers equation and KPZ-equation). In 
the limit of high Reynolds number, using the operator 
product expansion (OPE), Polyakov reduces the prob- 
lem of computation of correlation functions in the iner- 
tial subrange, to the solution of a certain partial differ- 
ential equation [19-20]. Yakhot recently [1,23] general- 
ize the Polyakov approach in three-dimensions and find 
a closed differential equation for the two-point generat- 
ing function of the "longitudinal" velocity difference in 
the strong turbulence (see also [21] about closed equa- 
tion for PDF of velocity difference for two and three- 
dimensional turbulence without pressure). On the other 
hands, recently [2] from detailed analysis of experimen- 
tal data of a turbulent free jet, R.Friedrich and J. Pienke 
have been able to obtain a phenomenological description 
of the statistical properties of a turbulent cascade us- 
ing a Fokker-Planck equation. In other words they have 
seen that the conditional probability density of velocity 
increments satisfy the Chapman-Kolmogorov equation. 
Mathematically this is a necessary condition for the ve- 
locity increments to be a Markovian process in terms of 
length scales. By fitting the observational data they have 
succeeded to find the different Kramers-Moyal(K.M) co- 



efficients and they find that the approximations of the 
third and fourth order coefficients tend to zero whereas 
the first and second coefficients have well denned limits. 
Then giving address to the implications dictated by [22] 
theorem they have got a Fokker-Planck evolution opera- 
tor. As an evolution equation for the probability density 
function of velocity increments, the Fokker-Planck equa- 
tion has been used to give the information on changing 
shape of the distribution as a function of the length scale. 
By this strategy the information on the observed inter- 
mittency of the turbulent cascade is verified. In their 
description and based on simplified assumptions on the 
drift and diffusion coefficients they have considered two 
possible scenarios in order to indicate that both the Kol- 
mogorov 41 and 62 scalings are recovered as possible be- 
haviors in their phenomenological theory. 
In this paper we derive the Kramers-Moyal equation 
from Navier-Stokes equation and show that how the 
higher order (n > 3) Kramers-Moyal coefficients tend 
to zero in the high Reynolds number limit. Therefore 
we find the Fokker-Planck equation from first principles. 
We show that the breakdown of the Galilean invariance is 
responsible for scale dependence of the Kramers-Moyal 
coefficients. Finally using the path-integral expression 
for the PDF we show that how small scale statistics af- 
fected by PDF's in the large scale and confirm the Lan- 
dau's remark that the large-scale fluctuations of turbu- 
lence production in the integral range can invalidate the 
Kolmogorov theory [12-13]. 

Our starting point is the Navier-Stokes equations: 



v t + (v • V)v = *A7 2 v *L + f( x ,t), 



(1) 



for the Eulerian velocity v(x, t) and the pressure p with 
viscosity v, in N-dimensions. The force f (x, t) is the ex- 
ternal stirring force, which injects energy into the system 
on a length scale L. More specifically one can take, for 
instance a Gaussian distributed random force, which is 
identified by its two moments: 

</„(*, *)/./(*.*') >=A(0)«(i-i')fc |u/ (x-x) (2) 

and < / M (x, t) >= 0, where /z, v = x\, X2, ••• , xn- The 
correlation function fc M , y (r) is normalized to unity at the 
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origin and decays rapidly enough where r becomes larger 
or equal to integral scale L. 

The force free N-S equation is invariant under space- 
time translation, parity and scaling transformation. Also 
it is invariant under Galilean transformation, x — > x + Vt 
and v — > v + V, where V is the constant velocity of the 
moving frame. Both boundary conditions and forcing can 
violate some or all of symmetries of force free N-S equa- 
tion. However it is, usually assumed that in the high 
Reynolds number flow all symmetries of the N-S equa- 
tion are restored in the limit r — > and r » 77, where n 
is the dissipation scale where the viscous effects become 
important. This means that in this limit the root-mean 
square velocity fluctuations u rms — V< v 2 > which is 
not invariant under the constant shift V, cannot enter 
the relations describing moments of velocity difference. 
Therefore the effective equations for the inertial-range 
velocity correlation functions must have the symmetries 
of the original N-S equations. For many years this as- 
sumption was the basis of turbulence theories. But based 
on the recent understanding of turbulence, some of the 
constraints on the allowed turbulence theories can be re- 
laxed [1]. Polyakov's theory of the large-scale random 
force driven Burgers turbulence [5] was based on the as- 
sumption that weak small - scale velocity difference fluc- 
tuations (i.e. \v(x + r) — v(x)\ « u rms and r « L), 
where L is the integral scale of system, obey G-invariant 
dynamic equation, meaning that the integral scale and 
the single-point u rms induced by random-forcing cannot 
enter the resulting expression for the probability density. 
According to [1] it has been shown that how the u rms 
enters the equation for the PDF and therefore breaks the 
G-invariance in the limited Polyakov's sense. We are in- 
terested in the scaling of the longitudinal structure func- 
tion S q =< (u(x + r) — u{x)) q >=< U q >, where u(x) is 
the x-component of the three-dimensional velocity field 
and r is the displacement in the direction of the x-axis 
and the probability density P(U, r) for homogeneous and 
isotropic turbulence. Let us define the generating func- 
tion Z for longitudinal structure function Z =< e xu >. 
According to [1] in the spherical coordinates the advec- 

tive term in eq.(l) involve the terms 0(3x57) , O(j^), 
0(j§p), O(-j^) [21]. It is noted that the advection con- 
tributions are accurately accounted for in equation of Z, 
but it is not closed due to the dissipation and pressure 
terms. Using the Polyakov's OPE approach Yakhot has 
shown the the dissipation term can be treated easily while 
the pressure term has an additional difficulty The pres- 
sure contribution leads to effective energy redistribution 
between components of velocity field and has non-trivial 
effect in the dynamics of N-S equation. Proceeding to 
find a closed equation for the generating function of Lon- 
gitudinal velocity difference, Z, the dissipation and pres- 
sure terms in eq.(l) give contributions and the longitu- 
dinal part of the dissipation term renormalizes the coef- 
ficient in front of the O(j) in equation for Z, [1]. Also it 
generates a term with order of 0(u) which can be writ- 



ten in terms of Z as A^j. Taking into account all the 
possible terms and using the symmetry of the PDF i.c 
P(U,r) = P(—U,—r), the following closed equation for 
Z can be found, [1], 



d 2 Z _ BodZ_ _ AdZ_ _ dZ 
dXdr X dr r dX dX 



3r 2 X 2 Z (3) 



where the parameter A, B and C to be determined from 
the theory. Also we suppose that has the struc- 

ture fc„ v ( r u) = fc (°)[ 1 - TF#<W - (riJ Y/' j) " 1 with 



fc(0) = 1 and = Xi — xj. The Gaussian assumption 
for " single — point" probability density fixes the value 
of coefficient C — and the C-term corresponds to 

the breakdown of G-invariance in the limited Polyakov's 
sense [5]. The A-term is responsible for interaction of 
transverse components of velocity field with the longi- 
tudinal component and produce an effective source and 
friction for the longitudinal correlation. 

In the limit r — > the equation for the probability 
density is derived from eq.(3) as, 



d_ v dP__ B dP _ _Ad_ 
dU dr dr r dU 



UP 



d 2 



L dU 2 



UP 



(4) 



Using the exact results S3 — — |er in the small 
scale, ( e is the mean energy dissipation rate) one finds 
A = 2±£, where B = -B > [1]. It is easy to 
see that the eq.(4) can be written as d r P = (—d{jU — 
B )- 1 [-(A/r)d u U+ (u rms /L)d1jU]P and so its solution 
can obviously be written as a scalar-ordered exponential 

/ dr' Lkm{U,t ') 

[22], P{U, r) - T(e+ ro P{U, r )) , where L KM 

can be obtained formally by computing the inverse op- 
erator. Using the properties of scalar-ordered exponen- 
tials the conditional probability density will satisfy the 
Chapman-Kolmogorov equation. Equivalcntly we derive 
that the probability density and as a result the condi- 
tional probability density of velocity increments satisfy a 
K.M. evolution equation: 



r) P r) n 

1 EH)"^P W M)P) 

n=l 



dr 



(5) 



Where D^ n \r,U) = ^U n + f3 n U n -\ We have found 
that the coefficients a n and (3 n depend on A and B, u rms 
and inertial length scale L which are given by the recur- 



sion relations We scale the velocities as U = 



and 



introduce a logarithmic length scale A = ln(^) which 
varies from zero to infinity as r decreases from L to n. 
Thus the form of DW (U, r) and £>( 2 ) (U, r) in the equiv- 
alent description would be D^(U,r) = —(tTb)U and 



DW{U,r) = {, 



-{2+B)(l + B))U (£) ' "™sV( 2+ _B), 

The drift and diffusion coefficients for various scales 
A, determined in the theory of Yakhot, show the same 



(1) 
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functional form as the calculated coefficients from exper- 
imental data [2,3]. 

In comparison with the phcnomcnological theory of 
Friedrich and Pienke we are able to construct a K.M. 
equation for velocity increments that is analytically de- 
rived from Yakhot theory which is based on just general 
underlying symmetries and OPE conjecture. Further- 
more this viewpoint on the equation (4) gives the ex- 
pressions for scale dependence of the coefficients in the 
K.M. equation. The important result is that scale de- 
pendent K.M. coefficients are proportional to u rms which 
gives a probable relationship between breakdown of G- 
invariance and scale dependence of the K.M. coefficients 
in the equivalent theory. The two unknown parame- 
ters A and B in the theory is reduced to one by fit- 
ting the £3 = 1, so all the scaling exponents and D^'s 
are described by one parameter, B. Considering the re- 
sults in [1,2] on which the value of B arc obtained, we 
have used the value B = 20 and have calculated the 
numerical values of KM coefficients. Ratios of the first 
three coefficients a n and /3„ are 0:3/012 = 0.04,04/02 = 

0.001, /3 3 //3 2 = 0.04, /V/3 2 = °- 001 - From the compari- 
son of numerical values of higher order coefficients we find 
that the series can be cut safely after the second term and 
a good approximation for evolution operator of velocity 
increments is a Fokker-Planck operator. According to 
[1] the value of parameter B = 20 is calculated numer- 
ically in the limit of infinite Reynolds numbers. Using 
this value for calculation of the numerical values of D\ 
and D 2 we find that the contribution of scale dependent 
terms are essentially negligible. This correspondence was 
referred to in [1] about the behavior of diffusion coeffi- 
cients for the limit of infinite Reynolds numbers. As it is 
well-known the Fokker-Planck description of Probability 
measure is equivalent with the Langcvin description writ- 



ten as [22], f£ = D^(U,X) + ^jD( 2 )(U,X)r](X), , where 
77(A) is a white noise and the diffusion term acts as a 
multiplicative noise. By considering the Ito prescription 
and using the path-integral representation of the Fokker- 
Planck equation we can give an expression for all the 
possible paths in the configuration space of velocity dif- 
ferences and thus demonstrate the change of the measure 
under the change of scale, i.e. 



P(U 2 ,X 2 \U ll X 1 ) = J V[U]e 



f A2 



4D 2 (U ,X) 



(6) 



When calculating, the measure of path integral is mean- 
ingful when some form of discretization is chosen [22] , but 
we have written it in a formal way. Using the forms of 
D 1 and D 2 and approximating them with scale indepen- 
dent ones in the infinite Reynolds number limit, one can 
easily see that the transition functional can be written in 
terms of In U. It is an easy way to see how the large scale 
A — > Gaussian probability density can change its shape 
when going to small scales A — > 00 and consequently give 
rise to intermittent behavior. Instead of working with the 



probability Functional of velocity increments the formal 
solution of Fokker-Planck equation as a scalar-ordered 
exponential [25], can be converted to an integral repre- 
sentation for the probability measure of velocity incre- 
ments when the D 1 = -ai{X)U and D 2 = a 2 (X)U 2 ,i.c 



P(U,X) = 



o7o(>) 



/+00 
1 
-00 



><f)(Ue 



71(A)- 



*)ds (7) 



where, 7o (A) = /„ (-ai(A') + 2a 2 (X'))dX' and 7l (A) = 
J A (-oi(A') + 3o 2 (A'))dA' and 7(A) = J* a 2 (X')dX' and 
4>{U) is the Probability measure in the integral length 
scales (A — > 0) . We consider the Gaussian distribution, 
(j>(U) = e~ mJj2 in the integral scale which is a reasonable 
choice ( experimental data shows that up to third mo- 
ments the PDF in the integral scale are consistent with 
Gaussian distribution [1]) and derive the dependence of 
the variance of the probability density on the scales in the 
limit when the original distribution satisfies the condition 
m <C 1. The result shows an exponential dependence like 
to — ► me 2 ^ where £ = 3o 2 — a\. The consistent picture 
with the shape change of probability measure under the 
scale is that when A grows, the width decreases and vice 
versa, which is reported in previous works as a simulation 
and experimental results. Moreover we should empha- 
size that the shape change is somehow complex which 
gives some corrections in order 0(m 2 [/ 4 ) even in this 
simplifying limit, i.e. m <C 1. Starting with a Gaussian 
measure at integral scales and using the calculated scale 
independent Fokker-Planck coefficients, we have numeri- 
cally calculated the PDF's for fully developed turbulence 
and Burgers turbulence in different length scales which 
their plots in Fig.[l] and Fig. [2] are completely compati- 
ble with experimental and simulation results [1,2,4]. The 
extreme case of Burgers problem (i.e.-B = 0) shows the 
ever localizing behavior as if in the limit of A — > 00 goes 
to a Dirac delta function which again is consistent with 
our knowledge about Burgers problem [1,4]. Clearly the 
eqs. (4) and (5) give the same result for multifractal 
exponent of structure function ,i.e. S n (r) = A n r^ n is de- 
rived to be £„ = [1]. 

In summary We have constructed a theoretical bridge 
between two recent theories involving the statistics of 
longitudinal velocity increment fluctuations in fully de- 
veloped Turbulence. On the basis of the recent theory 
proposed by V. Yakhot we showed that the probabil- 
ity density of longitudinal velocity components satisfy 
a Kramers- Moyal equation which encodes the Markovian 
property of these fluctuations in a necessary way. We are 
able to give the exact form of Kramers-Moyal coefficients 
in terms of a basic parameter in Yakhot theory B. The 
qualitative behavior of drift and diffusion terms are con- 
sistent with the experimental outcomes [2]. As the most 
prominent result of our work, we could find the form of 
path probability functional of the velocity increments in 
scale which naturally encodes the scale dependence of 
probability density. This gives a clear picture about the 
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intermittent nature in fully developed Turbulence. 

We should emphasize that the derivation of KM equa- 
tion is not restricted to the Polyakov's specific approach. 
One can show that similar results could be obtained by 
the conditional averaging methods [24-25]. Clearly ana- 
lytic form of the K.M. coefficients D^ n ' can be estimated 
numerically but analytic derivation is not possible [26]. 
Our work might be generalized to give a theoretical basis 
for the Markovian fluctuations of the moments of height 
difference in the surface growth problems like KPZ [18,27] 
and we believe that it would be possible to derive the 
Kramers-Moyal description for the statistics of energy 
dissipation [28]. 
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FIGURE CAPTIONS 



Figure 1. Schematic view of the logarithm of PDF in 
terms of different length scales. These graphs are numer- 
ically obtained from the integral representation of PDF 
at the Fokker-Planck approximation. The curves corre- 
spond with the scales L/r = 1.5, 2, 5, 10, 20. 

Figure 2. Schematic view of the logarithm of PDF 
in the Burgers turbulence (B = 0), in terms of different 
length scales. These graphs are numerically obtained from 
the integral representation of PDF at the Fokker-Planck 
approximation. The scales are L/r = 1.5, 2, 5, 10, 20. 
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